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Probing the Planck Scale with Proton Decay 
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We advocate the idea that proton decay may probe physics at the Planck scale instead of the GUT 
scale. This is possible because supersymmetric theories have dimension-5 operators that can induce 
proton decay at dangerous rates, even with i?-parity conservation. These operators are expected 
to be suppressed by the same physics that explains the fermion masses and mixings. We present 
a thorough analysis of nucleon partial lifetimes in models with a string-inspired anomalous U{l)x 
family symmetry which is responsible for the fermionic mass spectrum as well as forbidding _R-parity 
violating interactions. Protons and neutrons can decay via 7i-parity conserving non-renormalizable 
superpotential terms that are suppressed by the Planck scale and powers of the Cabibbo angle. 
Many of the models naturally lead to nucleon decay near present limits without any reference to 
grand unification. 



I. INTRODUCTION 

Baryon number is an accidental symmetry of the Stan- 
dard Model (SM). However, it is unlikely that baryon 
number will remain a symmetry up to the highest energy 
scales. Because of the gauge hierarchy problem it is 
strongly believed that the SM must be augmented by 
new physics at the TeV scale. There is no theoretical 
reason to believe that such new physics will still conserve 
U{1)b- In fact, baryon number violation is one of the 
necessary ingredients in models that explain the matter- 
antimatter asymmetry in the universe g. We therefore 
expect U{1)b to be violated at some higher energy scale. 
Current and future experiments of nucleon decay may 
thus be probes of high-scale physics. For example, 
the partial lifetime of the decay mode p — > K^v is 
greater than 1.6 x 10^^ years which places stringent 
constraints on scenarios of new physics. 

On the other hand, supcrsymmetry is considered the 
leading candidate for physics beyond the SM because it 
solves the gauge hierarchy problem once and for all and it 
is also consistent with grand unified theories (GUTs). It 
is often said that since quarks and leptons lie in the same 
GUT multiplets, an observation of proton decay will be 
a signal of a GUT. Indeed, the minimal supersymmetric 
SU (5) GUT has been excluded because the unification 



of gauge couplings forces the model into a region of 
parameter space where the proton decays too quickly . 

However, it is an under-publicized fact that super- 
symmetric models predict proton decay even without 
unification. Conventionally i?-parity (i?p) is imposed 
as an exact symmetry in order to prevent the proton 
from decaying through rcnormalizable operators. But in 
a generic supersymmetric model one can still write down 
i?p-conserving, yet baryon and Icpton number violating, 
Z? = 5 operators suppressed by a single power of the 
cutoff scale, which we take to be the reduced Planck 
scale Mpi ~ 2.4 x 10^® GeV. These operators come from 
the superpotential 
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Here Q,L are quark and lepton doublet superficlds 
and U,D,E are SU{2)l singlets. The SU{2)l indices 
are contracted between the terms in parentheses, while 
i,j,k,l are generational indices and a,b,c are color 
indices. From an effective field theoretic point of view 
one expects these operators to appear with coefficients 
Cl and Cr of 0(1). 

Supersymmetry's dirty little secret is that if one were 
to pick generic 0(1) numbers for the coefficients Cl and 
Cfl, the proton lifetime would be about 10^^ years, many 
orders of magnitude below the experimental limit. Thus 
the coefficients Cl and Cr must be highly suppressed in 
any viable supersymmetric model. 

We do not consider this embarrassment a death-blow 
for supersymmetric models, however. This is because 
the SM already contains highly suppressed dimensionless 



numbers, namely the Yukawa coupling constants. It is 
unsatisfactory to have an effective theory that is valid up 
to the Planck scale without explaining the origin of small 
coefficients for HuQiU j, H-pQiDj and H-pLiEj. Such an 
explanation should naturally determine the coefficients 
of QiQjQkLi and UiU jDkEi, likely suppressing them 
as well 1^]. In the context of supcrsymmetric models 
the most promising scenario for generating small Yukawa 
couplings is that of Froggatt and Nielsen |^ (see |^ for an 
exhaustive list of references), for which the suppression 
of proton decay has been demonstrated [Q, 0|- In ||] 
the tight experimental bounds on exotic processes were 
brought to bear on certain Froggatt-Nielsen models found 
in the literature. Indeed, many models were found 
to be incompatible with data solely due to insufficient 
suppression of QiQjQkLi- The relation between the 
fermion mass hierarchy in the SM and the suppression 
of proton decay also occurs in other non-supersymmetric 
models such as |§, ^ . 

In order to predict proton decay rates originating from 
Planck scale D = 5 operators we have to specify the 
framework of fermion masses and mixings. We would like 
the framework to be specific enough so that we can make 
quantitative predictions, but also general enough so that 
we can study the generic features of Planck scale proton 
decay. In this paper we will focus on a recently proposed 
class of models based on a single, anomalous Froggatt- 
Nielsen flavor symmetry jl^ . In this class of models the 
MSSM supcrficlds are charged under a horizontal U{l)x 
symmetry that is spontaneously broken by the non-zero 
VEV of a flavon field, A. The MSSM Yukawa terms are 
then suppressed by the ratio e = {A)/Mpi raised to the 
the appropriate power necessary to conserve U(l)x- The 
models in are ambitious: 

1. The U{l)x is the only symmetry beyond the SM 
gauge groups. 

2. The only fields charged under the SM gauge groups 
are those in the MSSM. However, there may be 
(hidden sector) fields charged only under U{l)x- 

3. The U{l)x charge assignments and breaking scale 
are inspired directly from string theory. Anomalies 
are canceled by the Green-Schwarz mechanism 
which places sum rules on the X-charges. 

4. The Cabibbo angle is calculated. It is set by 
the flavon VEV which is is determined by string 
theory, yielding the phcnomcnologically interesting 
value e = (A)/Mpi = 0.17-0.22. 

5. The charge assignments are chosen to yield the 
measured quark and Icpton masses, including CKM 
mixing. 

6. Neutrino mixings are also a result of the charge 
assignments. 

7. i?-parity is an exact (accidental) symmetry of the 
U{l)x charge assignments, preventing tree-level 
nucleon decay. 



That these seven goals can be achieved within a single, 
simple framework is non-trivial and encouraging. As 
was shown in [ pO[ , there is only a finite number of X- 
charge assignments with these properties.* It is very 
interesting that the predicted nucleon partial lifetimes 
in these models are generically within a few orders of 
magnitude of the current limits of about 10^^ years. This 
already rules out some specific charge assignments, and 
means that others will be directly tested very soon. It 
is important to note that in general these models do not 
have GUT-compatible X-charges, but they nonetheless 
lead to very interesting nucleon decay rates not much 
smaller than those predicted by GUTs. 

In this paper we systematically explore the nucleon 
partial lifetimes predicted by the U{l)x Froggatt-Nielsen 
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review these models, then in Section III 
quick estimates of the lifetimes. In Section lIVl we 
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describe our method of computing the nucleon decay 
rates. Section ^ contains our main results, namely that 
the U{l)x models of flavor are already constrained by 
current nucleon decay data, and they will continue to 
be tested by future experiments. Our conclusions are 
presented in Section while Appendix ^ carefully 
explains the dressing diagrams relevant for our analysis, 
and Appendix ^ discusses the effects of canonicalizing the 
Kahler potential and transforming from the interaction 
basis into the mass basis. 



II. FROGGATT-NIELSEN FRAMEWORK 

We will analyze the class of Froggatt-Nielsen models 
presented in [ p^ . In this section we summarize the 
results that are relevant for our analysis and refer readers 
to that paper for details. Below the reduced Planck 
scale the gauge group and particle content are those 
of the MSSM, with the addition of two right-handed 
neutrinos, a single flavon superfield A, and a generation- 
dependent U{l)x gauge group with anomalies canceled 
by the Green-Schwarz mechanism fll|| . Since all matter 
fields are charged under U{l)x, their couplings to the 
flavon are determined by U{l)x invariance. The goal 
is to have the flavon couplings generate the generation 
dependent masses and mixings in the fermion spectrum. 

For example, masses for the up-type quarks originate 
from the non-renormalizablc operator 



gu 



A 



PI 



Xh^^+Xq.+X- 



Hu Qi U j 



(2) 



where the flavon charge is defined to be Xa = — 1 and 



With further assumptions, discussed in jlo| , these models can 
also explain the scale of neutrino masses by invoking only two 
right-handed neutrinos and using only two fundamental mass 
scales, A/pi and rrigoft ~ 1 TeV. 
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gu'^^ are 0(1) cocfBcicnts that are undetermined within 
the effective theory below Mp\. The dimensionless gu^^ 
is zero if Xhu +XQi +^Fj negative or fractional. The 
powers of A in Eqn. ^ compensate the U{\)x charges 
of Hu Qi Uj yielding U{l)x gauge invariants. The Dine- 
Seiberg-Wen-Witten-generated Fayet-Iliopoulos term for 
U{l)x induces a VEV for the scalar component of the 
fiavon field [l5| . The VEV is given by (A) = 

eMpi, where e naturally turns out to be the size of the 
sine of the Cabibbo angle, e ^ 0.2. After U (l)x-breaking 
one gets the Yukawa couplings 



Gn 
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(3) 



which are suppressed by powers of e. In this way the 
charge assignments of the MSSM matter fields determine 
the fermion mass hierarchy. Once the X-charges arc cho- 
sen to reproduce the mass spectrum, the e-suppressions 
of other superpotential terms are determined up to 
0(1) coefficients, with higher-dimensional operators like 
QiQjQkLi being suppressed by additional powers of Mpi. 
It is this connection that we will be exploiting to estimate 
the nucleon decay lifetimes within this Froggatt-Nielsen 
framework. 

The X-charges are determined by requiring cancel- 
lation of chiral anomalies between U{l)x and the SM 
gauge group, prohibition of all i?p- violating interactions, ^ 
and generation of the phenomenologically viable fermion 
masses and the CKM matrix. The X-charges of the 
MSSM superfields are then specified by six integer pa- 
rameters, as demonstrated in [ p^ and displayed in their 
Table 2. These X-charges are further constrained by the 
phenomenology of neutrino mixing, which fixes two of 
those six parameters . The X-charge of each superfield is 
determined by the four remaining parameters, x, y, z and 
A^, and is shown in Table |. The physical significance 
of these parameters will be explained below. 

Only the three parameters x, y and z are relevant 
for nucleon decay. They are restricted to a small set of 
integers, which leads to 24 distinct models with different 
nucleon decay signatures. The allowed values are x = 
0, 1, 2, 3; 2/ = -1, 0, 1, * and z = 0, 1. 

The parameter x is related to the ratio of the bottom 
and top quark masses, and is thus also connected to tan /3, 
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i + Xh^ + A" 
Xl^ — Xli — z 
-Xhj, + 4 ~ Xli + X + z 
~Xht, + 2-Xi,i +X + Z 
—Xh-o — Xli + X + z 



TABLE I: i?p-conserving A-charge assignments of the MSSM 
superfields in terms of four parameters: a;^, z, . The 
charge assignments are as in Table 2 of lid] with their 
and A31 both taken to be —z as required by neutrino mixing 
phenomenology. 



Since e ~ 0.2, larger x corresponds to smaller tan/3, 
but because of unknown 0(1) coefficients, we cannot 
determine tan /3 exactly. To simplify our analysis we 
will choose a specific, reasonable value of tan /3 for each 
value of X, namely tan/3 = 50, 20, 5, 3 for x = 0, 1, 2, 3, 
respectively. 

The three choices for the parameter y determine the 
texture of the up- and down-quark Yukawa matrices and 
therefore also the CKM matrix. One finds 
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(4) 



y = Q 



VcKM 



(6) 



t All Rp-even terms are supposed to have an overall integer X- 
charge, while all ijp-odd terms are supposed to have an overall 
half-odd-integcr X-charge. 

t In Ref. y = —7, —6 are also considered because at first sight 
they give a viable CKM matrix once the supersymmetric zeros 
are filled in. However, Ref. [h6| demonstrates that this is not 
the case. See also Refs. [h7|. Tncrcfore we do not consider y = 
—7, —6 in this paper. 
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Finally, the parameter z is related to the charged 
lepton mass spectrum. 



m„ 



(8) 
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In combination with neutrino phenomenology and the 
requirement that Rp is conserved by virtue of the X- 
charges, it turns out that z specifies the texture of the 
MNS matrix: 



Ui 



'MNS 




(9) 



Here we see that z ~ Q corresponds to an "anarchical" 
MNS matrix [l^, whereas z = 1 corresponds to a 
"semi-anarchical" MNS matrix [ pO[ . 

The minimalist approach is to add only two right- 
handed neutrinos. When this approach is taken the 
fourth parameter, A^, can take two different values for 
each set of {x^y,z}, but it has no impact on nucleon 
lifetimes and so it is irrelevant for our purposes. However, 
our analysis does not depend on the number of right 
handed neutrinos. 

In addition to specifying the CKM and MNS textures, 
the parameters {a;,y,z} also have bearing on other as- 
pects of the UV physics. For example, SU{5) invariance 
is only consistent with y = 1 and z = 0. Also, z ~ 
1 prohibits the operator HxiHy, thus allowing the /.j- 
term to be generated by the Giudice-Masiero mechanism 
If one has a preference towards a specific 
physical scenario, the number of models to choose from 
is significantly reduced. For example, the models in ||lo[| 
were chosen to have z = 1 for the sake of a natural 
/i-parameter, ?/ = to obtain the most natural CKM 
matrix, and a; = 2, 3 to avoid proton decay limits. 

Other criteria may be used in choosing one's favorite 
model. For instance, the values of the X-charges are 
much more aesthetically pleasing in some models than 
others. For example, taking x = 2, y = l, z = and 
one finds the X-charges of all MSSM superfields 
are integers or half-odd-integers, as shown in Table ||. 
However, in this paper we will treat all 24 distinct models 
equally and focus only on their predictions for nucleon 
decay. 
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-3/2 



X 



-3, X 



Hu 



TABLE II: An example of an X-charge assignment for the 
matter superfields in the model with x = 2, y = 1, z = Q, and 
= 9. helps to determine the specific charges, but its 
effects cancel out in all the nucleon decay operators. 



III. QUICK ESTIMATES 

A quick estimate shows that the operator Q1Q1Q2L3 
leads a proton lifetime of ^ lO^^e"^" years, where n = 
9 + X + y + z is the sum of the A"-charges for the four 
superfields. Since e ^ 0.2, we need to have n > 11 to 
satisfy the experimental limit of about lO'^^ years for p 
K^Vr- This demonstrates the general trend that larger 
values of x, y, and z will lead to greater suppression of 
the decay rate and hence a longer lifetime. 

It turns out that the nucleon lifetimes predicted by 
our Froggatt-Nielsen model are in the same ballpark as 
those predicted by GUT models. This can be understood 
relatively easily. For the operator Q2Q2Q1L2 we can 
rewrite the suppression in terms of Yukawa couplings 
(ignoring 0(1) factors). 



1 
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PI 



Mpi 1713/2 
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(10) 

where we have used the phenomenological constraints 
on the X-charges from Section 5 in |]lOj to identify the 
Yukawa couplings. This is to be compared with the 
coefficient of the same operator coming from an SU (5) 
GUT theory, which is 



1 



M, 



GUT 



22 22 \ 
- (-T/y Ue Ac . 



(11) 



Equating e with Ac and noting that the /x-parameter can- 
not be too different from the gravitino mass, the differ- 
ence between the coefficients is simply Xc {Mqut / Mp\) . 
Considering the fact that these two models involve 
completely different physics, we find this difference inter- 
estingly small. The estimation above is obviously very 
crude and calls for a more rigorous analysis which is the 
goal of this work. 



IV. NUCLEON DECAY COMPUTATIONS 

In this section we describe our calculations of the various 
nucleon decay modes. We use the standard methods 
for the computations, following Q and j23|. The 
starting point is the higher-dimensional UV Lagrangian 
shown in Eqn. (|l|). The dimension-five fermion-fermion- 
scalar-scalar terms arising from these operators can be 
"dressed" , whereby the scalars exchange a gaugino or 
higgsino and are converted into fermions, as shown in 
Figure |l|. Below the scale of the soft masses they become 
four-fermion operators and contribute to nucleon decay. 
The matrix elements of these operators can be evaluated 
using the well-known chiral Lagrangian techniqu e Jp4| . 

The quark doublet superfields Q in Eqn. (|^)are 
taken to be in the SuperCKM basis, namely Q = 
([/, D') = (U, VckmD) where U and D are the mass 
eigenstates. Thus the couplings to the mass eigenstate 
quark operators will come with various CKM factors. 
However, we cannot necessarily neglect operators with 
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FIG. 1: An example of a dimension-5 operator that can be 
dressed by neutral or charged winos and higgsinos. 



small CKM factors because they may be offset by the 
lesser degree of e-suppression in C]^"'^. In other words, 
operators QiQjQkLi containing third generation quarks 
will have small CKM couplings to the first generation 
quarks in the nucleon, but will have correspondingly 
larger coefficients C^^^^ due to the X-charge assignments 
for the third generation. 

Note that the U{\)x basis is not necessarily the same 
as the SupcrCKM basis nor the mass eigenbasis. Thus it 
is important to take these various changes of basis into 
consideration when determining the e-suppression of the 
nucleon decay operators. For the most part this only 
effects the 0(1) coefficients, which we are ignorant of 
anyway. Nevertheless, we have performed a thorough 
analysis of this effect as described in Appendix 

The operators QiQjQkLi and UiUjDkEi are not all 
independent. Fierz identities may be used to convert 
all of the operators into a smaller set which we choose 
as our basis of independent operators. One can further 
reduce the number of operators in this set by noting 
that contributions to nucleon decay can only come from 
operators with at least one first-generation superfield. 
The Fierz identities and a list of independent operators 
are given in Appendix Since we do not know the exact 
coefficients in front of these operators, wc will assume 
that in the UV all of the operators in our basis are 
generated and that their unknown coefficients are given 
by the e-suppression determined by the AT-charges and 
0(1) numerical factors. 

As argued above, the e-suppression is essential for the 
predicted nucleon lifetimes to be above the current limits. 
This suppression is determined by the X-charges of Qi, 
Ui, Di, Li and Ei, which are in turn determined by the 
choice of the model parameters x, y, and z. Note that e 
itself also has mild dependence on the model parameters. 



.93 



47rV2 



V3(6-f x + z) 



(12) 



where the QCD gauge coupling 53 is roughly 0.72 at 
high energies. The dependence of the e-suppression on 
the choice of X-charges leads to the different patterns 
of nucleon partial lifetimes that we study in the next 
section. 



Rcnormalization group effects enhance the nucleon 
decay operators due to running between Mpi and rup. 
For simplicity, we compute the supersymmetric running 
from Mpi down to mz from the gauge couplings, ignoring 
corrections proportional to Yukawa couplings. The left- 
and right-handed operators, QiQjQkLi and U iU jD^Ei, 
are enhanced by a factor 



Here 



ai(Mpi) 
aiimz) 



Oi2{mz) J \as,{mz) 



(13) 



^[mz) - ^\n{q/mz) with bi = 
■- (1/33,3,-4/3) for LH and 



(33/5,1,-3) and c, 

(2/11,0,-4/3) for RH operators. Also, gf = (5/3)g'2 
Numerically this enhances the LH operators by a factor 
of 9.3 and the RH operators by a factor of 5.6. 

In general there are many different ways the dimension- 
five operators can be dressed to give four-fermion oper- 
ators. In our scenario we assume that supersymmetry 
breaking is sufficiently flavor blind so that the flavor mix- 
ing in giuino-quark-squark vertices is negligible. Then in 
the limit of degenerate squarks and degenerate slcptons 
the contributions from gluino dressing cancel due to a 
Fierz identity, as explained in more detail in Appendix 
There it is also shown that bino dressing does not con- 
tribute to nucleon decay. All contributions from neutral- 
higgsino dressing are extremely suppressed by small, 
first-generation Yukawa couplings, so we ignore them. 
We do include the leading contributions from charged- 
higgsino dressing. In the case of chargcd-lcpton decay 
modes these contributions turn out to be negligible, 
but p3| has shown that they can be important for the 
antincutrino decay modes. The bulk of the contributions 
come from charged-wino dressing, though the neutral- 
wino dressing can be comparable for the antineutrino 
decay modes. 

The dressing leads to a four-fermion vertex which must 
then be run from the scale of the soft masses, roughly 
10^ GeV, down to 1 GeV where they mediate nucleon 
decay. For simplicity wc compute the QCD running only 
from mz down to 1 GeV, which leads to an enhancement 
given in |2^ ], namely 



0^3 (Mow) 

a;3(A^high) 



2/Po 



(14) 



where /3o — 11 — |/ and / is the number of quark flavors 
with masses less than Afhigh- There are three of these 
enhancement factors for the successive steps between 
mz, JTifc, mc, and mp. In each step the QCD coupling 
is given by as (g) — 47r/(/3o ln((7^/A/^)), where Af is 
determined by requiring that as be continuous between 
successive intervals. Numerically we use mf, = 4.2 GeV, 
TOc = 1.2 GeV, and mp ~ 1 GeV, which gives a total 
contribution of 1.3, and is quite insensitive to the exact 
thresholds for mi, and m^. Together with the high scale 
running mentioned above, there is a net enhancement of 
approximately 12 for Cl and 7 for Cr. 
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The chiral Lagrangian can be used to determine the 
decay rate induced by the four-fermion operators [^ . 
The partial decay width of a nucleon Bi into a meson 
Mj and lepton Ik is given by§ 

rw^M,o^3^(i-5;y(i.fr+i.4fr) 

(15) 

where A^*^ can be found in Table 1 of [p3| , given in terms 
of the coefficients of the various four-fermion operators 
generated by dressing of the superpotential operators 
shown in Eqn. (|l|). The numerical values of A^^ft depend 
on several input parameters. The coupling between 
baryons and mesons in the chiral Lagrangian are given 
by D and F. As in we take D = 0.81, F = 0.44, and 
,U = 139 MeV. 

Our computations of nucleon lifetimes are subject to 
two types of uncertainties, those that reflect our current 
lack of knowledge and those that are inherent to the 
effective theory framework of our model. The first type 
includes the uncertainty in the hadronic matrix clement 
and the masses of the superpartncrs, both of which 
we assume will be pinned down more precisely by the 
time the next generation nucleon decay experiments start 
taking data. The second type includes the unknown 
0(1) coefficients that come with the non-renormalizable 
superpotential operators, which cannot be determined 
without a full theory of the Planck scale physics. We 
will discuss each of these uncertainties in turn. 

The value of the hadronic matrix element must be eval- 
uated to convert the four-fermion operators into nucleon 
lifetimes. In the framework of the chiral Lagrangian 
this appears as two parameters /3p and ap, which are 
related by ap — ~f3p. Unfortunately, the value of Pp is 
only known to roughly a factor of 10. The conservative 
range often taken is Pp = 0.003-0.03 GeV^. Recently 
there has been some progress evaluating (3p on the lattice, 
yielding the results Pp = 0.014 ± 0.001 GeV^ in ||l and 
Pp = 0.007± 0.001 GeV^ in Note that the errors are 
only statistical, and do not reflect systematic effects such 
as quenching. Since Pp appears in the amplitude, any 
uncertainty in its value gets squared in the decay rate or 
lifetime. However, this uncertainty is common across all 
modes and all different models, so it simply represents an 
overall rcscaling of the lifetimes. For our computations 
we will use an intermediate value, Pp = 0.01 GcV"^. A 
smaller value of Pp would give a smaller decay rate and 
a longer lifetime. Allowing for the range 0.003-0.03 thus 
corresponds to a change in our results by a factor of 10 in 
both directions. Again, we can look forward to reduction 
in the uncertainty in the hadronic matrix clement, which 
will reduce the uncertainty in our computations. 



^ Eqn. ( p^ is written in the notation of p^ : fc is a generational 
index, out j refers to tt'', r?", /C", tt^ or K'^ , and i refers to 
proton or neutron. 



The scale of the superpartner masses is another un- 
certainty in our computations that will be reduced if 
supersymmetry is discovered at a collider. The super- 
partner masses enter through the loop in the dressed 
diagrams such as the one in Figure |l|. We assume 
that the squarks and sleptons all have a common mass, 
TOsoft- If the gaugino (or higgsino) mass is much less 
than the squark mass, the contribution from the loop 
is given by Mgaugino/'T^softJ whereas if the gauginos are 
degenerate with the squarks at msoft , then the loop factor 
is l/(2msoft). One can imagine extremes where all the 
superpartncrs are relatively light, near 100 GeV, leading 
to the shortest lifetimes, or the opposite extreme where 
the scalars are heavy (TeV) and the gauginos are light 
(100 GeV), leading to the longest lifetimes. 

We can cover the different possibilities by taking all 
superpartncrs to be degenerate at rrisoft and then allowing 
TOsoft to range from 100 GeV to 10 TeV. The latter case 
docs not mean the superpartner masses are actually 10 
TcV, but rather represents the case where squark and 
slepton masses are around 1 TeV and wino masses are 100 
GeV. In our computations we will choose a middle ground 
by assuming of TOsoft = 1 TeV. With this assumption, 
the proton lifetime simply scales inversely with to^qj-^.. 
Choosing one of the other two scenarios described above 
would change the nucleon lifetimes by a factor of 100. 
We reiterate that this uncertainty will disappear once 
the superpartner spectrum is measured, most likely at 
the LHC. 

Our ignorance of the 0(1) coefficients appearing in 
front of each operator in Eqn. (|l|) is an uncertainty 
inherent in the framework of effective field theory. A^^^i 
in Eqn. ( p^ ) are each the sum of several contributing 
amplitudes from different operators that can interfere 
with one another. However, since we do not know either 
the exact magnitudes or phases of the 0(1) coefficients 
in the UV Lagrangian, we cannot know whether the 
different contributions will interfere constructively or 
destructively. As an extreme example, one could pick 
a fine-tuned set of numbers such that the amplitude for 
a certain decay mode vanishes completely. However, this 
is unlikely to occur in the real world unless there is some 
symmetry of the high-energy theory that enforces such a 
cancellation. 

To remove the effects of such unlikely cancellations 
from our computations we will add the various terms 
contributing to Al^r in quadrature. This effectively gives 
the average of many computations with different random 
phases between the individual contributions. We take 
this as the central value of our results. However, it is 
possible that even without fine-tuning there could be 
significant effects due to interference between amplitudes. 
In order to estimate this effect we perform the calculation 
in two other ways, plotting the results as upper and 
lower error bars around the central value. The lower 
error bar is determined by choosing all the amplitudes to 
have the same phase, thereby interfering constructively. 
The upper error bar is determined by choosing the 0(1) 
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numbers to all be +1 or —1 such that amplitudes of 
similar sizes interfere destructively. Even though these 
are not rigid upper and lower bounds, they give a 
good sense of the possible variations due to interference 
between amplitudes. However, it is important to note 
that these are in no way Icr error bars. 

Finally, we will compare all of our computations of 
the partial lifetimes to the experimental limits, which 
are taken from 101 and |28l and are shown in Table [II. 



Mode 


Expt. Limit 


p — > i\ iy 


ID 


p — > 7r''"'/i^ 


o / 


p ^ n. ^ 


1 n 


p 


7.8 


p ^ 7v°e+ 


50 


p ■k'^v 


0.25 


p ^ K°e+ 


5.4 


p —* rf'e^ 


11 


n K°iy 


3.0 


n — + rPv 


5.6 


n — + -K^v 


1.1 


n —* Tv^ fi'^ 


1.0 


n TT^e^ 


1.6 



TABLE III: The various nucleon decay modes we consider, 
together with the experimental upper bound on the partial 
lifetime. Experimental values are from Q| and [ psf . 



V. RESULTS 

In this section we present the results of our analysis of the 
nucleon partial lifetimes within the context of the U{l)x 
flavor model of When the various amplitudes are 

added incoherently, as described in the previous section, 
the most constraining mode for all models is p — > 
because it has both a relatively large rate and a stringent 
experimental bound. First we will focus on this mode and 
then come back to the more promising models and look 
at other decay modes in more detail to see what we can 
expect to learn from nucleon decay experiments in the 
near future. 

As described in Appendix the Froggatt-Nielsen 
suppression of the operators QiQjQkLi and U iU jD^Ei 
are not trivially determined by U {l)x conservation, since 
they are affected by the canonicalization procedure of 
the Kahler potential and the transformation from the 
interaction basis into the mass basis. In Appendix |^ we 
will show that this is not an issue. 
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x=Q x=\ x=2 x=3 x=0 x=l x=2 x=3 
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FIG. 2: Plot of proton partial lifetime in years for the mode 
p — > K^u, shown for models with y = —1, 0, 1. Models with 
z = Q (1) are shown on the left (right), and within each half 
X increases from left to right, corresponding to a decrease in 
tan/3. The error bars are not la bars. Rather, the central 
point corresponds to incoherent addition of amplitudes. The 
upper and lower ends of the 'error bar' are the lifetimes 
computed by adding the various amplitudes destructively and 
constructively, respectively. The horizontal line shows the 
experimental lower limit of 1.6 x 10"^"^ years. The calculation 
was done with m-soft = 1 TeV and flp = 0.01, but the scales 
on the right show the overall shift that would result from 
changing either mgoft or f3p. 



In Figure ^ we show the partial lifetimes in the p 
K~^D mode for the 24 models, using the three different 
methods of adding the different contributions: incoherent 
(central value) , constructive and destructive (tips of error 
bars). The plot is divided in half, with z = on the left 
half and z = 1 on the right half, and in each half the 
value of X increases from left to right, corresponding to 
a decrease in tan [3. The current experimental bound is 
plotted as a horizontal line. We have included scale bars 
at the right to indicate the size of the overall systematic 
uncertainty due to mgoft and Pp. 

The first thing to notice is that many of the models 
are disfavored, unless they arc fine-tuned to give a highly 
destructive interference or the supersymmctric spectrum 
is carefully chosen. The models that survive best are 
those with a high x which corresponds to a lower tan /?. 
In fact, of the models with a; = or a; = 1, only one 
model has a prediction above the current experimental 
limit even when the amplitudes are added destructively. 
Recall, however, that the uncertainties in f3p and mgoft 
are not included in the error bar, which can potentially 
change the prediction by the factors shown graphically 
to the right of Figure g. Nevertheless, it seems clear 
that models with a; = or a; = 1 are disfavored by 
current proton decay limits, barring extremely delicate 
accidental cancellations, as was anticipated in based 
on a rough estimate. For this reason we will focus on 
models with a; = 2 or a; = 3 in what follows. From 
Figure H we also see that a long proton lifetime favors 
a semi-anarchical MNS matrix {z — I) and not one 
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determined by anarchy (z = 0). Note that the preference 
for z = 1 is also consistent with the Giudice-Masiero 
mechanism that naturally produces a fi-teim of the right 
size. However, this preference is only mild (within the 
uncertainty) . 

The current and upcoming experiments are expected 
to be able to detect charged leptons with a high efhciency. 
It is thus interesting to also focus on modes with charged 
leptons in the final state. In fact, for the proton the next 
modes to appear after p K^D arc generally p ir'^e'^ , 
p 7r°/^+ , and p K'^fi^ . In Figure || we specialize to 
models with x = 2 or 3 (i. e. tan p < 10) and show the 
expected lifetime from all four of the above decay modes. 
Here we plot the lifetime computed using incoherent 
addition of amplitudes, but the uncertainty due to 
unknown 0(1) numbers should be kept in mind. We 
see that most models which survive the constraint from 
p K'^v have a lifetime for p ii^ ji^ that is within 
two or three orders of magnitude of the experimental 
bound, while p is only slightly larger, and 

p — > 7r°e+ can potentially be smaller. This raises the 
exciting possibility of two or three decay modes being 
detected in the coming round of experiments. Figure ^ 
also shows how the relative decay rates can be used to 
distinguish models with different phenomenology. For 
instance, all models with z = (anarchy) have p —> n'^e'^ 
as the dominant decay mode with a charged lepton in the 
final state, whereas its lifetime is substantially higher for 
the models with z ~ 1. In both cases p tt'^h^ and 
p —> if ^/i"*" are comparable. 

To further illustrate the differences in the pattern of 
decay modes from different models, we will focus on three 
specific models, namely Model 1 with x = 3, y ^ —I, z = 
1, Model 2 with x = 2, y = 1, z = 0, and Model 3 
with X = 3, y ~ 0, 2 = 1. Model 3 was specifically 
studied in ||l^ . Model 2 is has a particularly nice charge 
assignment as shown in Table ||. In Figure ^ we show 
the normalized partial lifetimes for these three models 
in eight proton decay modes (left side) and five neutron 
decay modes (right side). 

Various modes, if discovered, may serve as a discrim- 
inator between models. For example, among the three 
models shown in Figure ^, any mode involving a muon 
in the final state can differentiate Model 1 from Models 
2 and 3. The decay mode p n^e^ can differentiate 
Model 2 from Model 3, etc. 



VI. CONCLUSION 
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FIG. 3: Comparison of proton lifetime in years for four 
different decay modes. The upper plot shows the computed 
lifetime for p — > K^v (x, left axis) and p K^'/i'^ (A, 
right axis). The lower plot shows p — » 7r"e^ (•, left axis) 
and p — > TT^jj.'^ (□, right axis). For each plot the data for 
the two modes have been scaled so that the experimental 
limits, represented by the horizontal line, coincide, hence the 
different scales on the left and right axes. The lifetimes are 
plotted for all twelve models with tan/3 < 10. 



generically right near the current experimental limits, 
even without Grand Unification. In fact, current bounds 
disfavor many of the 24 distinct models of this type. 

The main unknown in this program is whether super- 
symmetry exists, and if so, what the mass spectrum is. 
Data from the LHC and a Linear Collider, along with 
improved lattice calculations of the matrix element, will 
remove much of the uncertainty (see Fig. ^ leaving us 
with much tighter constraints. Furthermore, discovery 
of several proton decay modes would serve as a good 
discriminator between the various model parameters. 
Similarly, a measurement of ^13 in neutrino experiments, 
or of tan /? at a collider, can help narrowing down the 
choice of viable models. We thus conclude that upcoming 
experiments may truly be probing physics at the Planck 
scale. 



We have shown that the search for nucleoli decay is a 
powerful probe of physics at the Planck scale. We focused 
on a a class of ambitious, string motivated, Froggatt- 
Nielsen models. These models explain the masses and 
mixings of all SM fermions while automatically enforcing 
i?-parity as an accidental symmetry. In the context of 
these models, we have shown that operators suppressed 
by the Planck scale lead to nucleon decay rates that are 
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APPENDIX A: DRESSING DIAGRAMS 

The superpotential operators QiQjQkLi and UiUjDkEi 
are not all independent. In particular. 



QiQjQkLi — QjQiQkLi 
QiQjQiLi = --QrQiQjLi (no sum) 
Ujj.DkEi - -U,U,DkEi. 



(Al) 



Furthermore, contributions to nucleon decay can only 
come from operators with at least one first generation 
superfield. The list of independent operators that con- 
tribute to proton decay is shown in Table IV. Since 



we do not know the exact 0{1) coefficients in front of 
these operators, we will simply work with an indepen- 
dent basis of operators and assume that their unknown 
coefHcients are given by the e-suppression determined by 
the X-charges times 0(1) numerical factors, excluding 
accidental cancellations. The operators are then dressed 
by a gaugino or higgsino as shown in Figure |l|. In general, 
there are many ways each diagram can be dressed. 
However, it has been long known that in the limit of 
degenerate squarks many dressing combinations cancel 
due to a Fierz identity p9|, 30 . Here we demonstrate 



Operator 


Independent set of (ijkl) 


QiQjQkLi 


(1120, (1130, (2210, (3310, (1230, (1320, (2310 
(1210, (1220, (1230, (1310, (1320, (1330 



TABLE IV: Independent operators relevant for nucleon decay 
as labeled by their generation indices. 



this cancellation explicitly using gluino dressing as an 
example, following the description in pl| . Because 
gluinos are flavor blind, they couple to up- and down- 
type squarks equally and do not change flavor. Since 
they do not couple to leptons, there are three ways that 
gluinos can dress a given operator UiU jDkEi. 



/([/„[/,) {E^DkXmU^) 
f{U,,Dk) {EiUj){DkU,) 
f{Dk,U,) {EiUi){U,Dk) 



(A2) 



In these expressions the fields arc Weyl spinors con- 
tracted within the parentheses. The function /(x, y) 
comes from the loop integral and depends on the masses 
of the gluino and the squarks x and y. Thus f{x,y) 
encodes all the flavor dependence, so in the limit of 
degenerate squarks each factor of / above is common 
to all three terms. The interesting fact is that the sum 
of those three operators vanishes by the following Fierz 
identity: 

{AB){CD) + {AD){BC) + {AC){DB) = 0, (A3) 

which can be easily shown by rewriting the epsilon 
tensors used to contract the pairs of Weyl spinors in 
terms of Kronecker delta functions. Thus the sum of the 



three operators in Eqn. (A2) will vanish whenever they 
all have the same coefficient out front, which occurs in the 
degenerate squark limit. This can also be explained in 
words by noting that in the final four-fermion operator 
the three quark fields must be antisymmetric in color, 
antisymmetric in flavor (not generation, but flavor) and, 
since they are fermions, antisymmetric in spin as well 
to make the total operator antisymmetric under fermion 
exchange. But since there are three quark fields but only 
two spin states for each Weyl spinor, there is no way to 
make an operator that is completely antisymmetric in 
spin, hence the whole collection must vanish. 

This argument holds for all the gluino dressings of 
QiQjQkLi and UiUjDkEi. In addition, it also holds for 
bino dressing of QiQjQkLi | |3l[ |. The bino and neutral 
higgsino dressings of UiUjDkEi do not lead to any 
operators that contribute to nucleon decay, due to the 
inevitable presence of 2nd or 3rd generation up-quarks. 
Thus the only relevant dressing of the SU{2) singlet fields 
in UiUjDkEi is by charged higgsinos. The QiQjQkLi 
operators get contributions from charged and neutral 
winos and higgsinos. In general, the dominant contri- 
butions are usually from wino dressing of QiQjQkLi , 
except for the UiUjDkEi contribution to the Ki/t final 
state when tan/3 is large, as pointed out by [p3[. 
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APPENDIX B: 



TRANSFORMING INTO THE 
MASS BASIS 



Once the flavon acquires a VEV the coefficients C]f ^ are 
each determined by a dimcnsionlcss 0(1) coefficient times 
an e-suppression. In principle, some coefficients may be 
exactly zero due to negative overall X-charge, the so- 
called "supersymmetric zeros", but this does not occur 
in the models we consider. The naive C];'^ cannot be 
directly plugged into Eqn. (p^). First one needs to take 
into account two superfield-transformations, namely the 
canonicalization of the Kahler potential (CK), and then 
the transformation into the mass basis of the quarks and 
leptons (TM): 

1. From the outset the Kahler potential need not 
have the canonical form (sec g, ^ |4[ ^ for 
details). For example, the Kahler potential for the 
quark doublets takes the form 



Q, Q, = Cq.Q 



Cq-Q 



(Bl) 



C 



ijkl 



CK+TM 'J*-'' 



As we lack knowledge of the 0(1) coefficients, this 
can only be done approximately, supposing no accidental 
cancellations. One can show |34l that 



Cq 



-X 



Q3 



etc. 



(B5) 



from which it follows that the CK fills up supersymmetric 
zeros but does not change the e-suppression of any other 
nonzero entries. Using Table ^ this gives 



( 1 

g|l+!/l 

y g|3+al 



1 £2 



(B6) 




(B7) 




(B8) 



where Hq is a Hermitian matrix with hierarchical 
entries, generated when the flavon acquires a VEV. 
It can be diagonalized by the matrix Cq. This 
redefinition, i.e. the "CK", affects the coupling 
constants of the superpotential, e.g. 



Gi 



u 



Gi 



CK 



1 



U 



Cq-^'^ ■ Gu ■ (B2) 



2. Transforming the superfields to the mass basis, the 
coupling constants are then subject to the "TM" . 
One has 

qJ^k ^ Uul^ ■ diag{mu,mc,mt} -Uij^, (B3) 

etc., the 14,,, being unitary. 

The other superpotential coupling constants have to be 
transformed correspondingly, e.g. 



iU-r 



k r 



c 



'i' j' k' V 
R 



(B4) 



(B9) 
(BIO) 

The U,,, are given in Table 0, having used the expressions 




1 36 



Since the C]^*^ we consider do not have any 
supersymmetric zeros, their e-dependence is not changed 
by the CK. As for TM, comparing Table with the C... 
given above, one finds e.g. Uul * ^ Cq and Ujj^^ ^ C^j. 
Since Cq^ ^ Cq (and likewise for the others), we find 
that the C]^^ are not changed substantially when going 
into the mass basis either. Thus the naively calculated 
e-suppression for C];'^ remains the same after changing 
to the mass basis. 
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1 1 e'+y e'+y\ 

\^3+y ^2 ^ J 


( 1 e^"^ \ 
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/ 1 \ 

e"^ 1 1 
V 1 1 / 


/ 1 \ 
1 
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the five pairs of quark mass matrices displayed in The ones for the leptons (y-independent) are only valid for the X-charge 
assignment used in this article. Unj^ is the almost unitary matrix which Schur-diagonalizes the symmetric mass matrix of the 
light neutrinos. Note that these approximate TM matrices reproduce the CKM matrix and the MNS matrix, as they should. 
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